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Effect of chemistry and turbulence on NO 
formation in oxygen- natural gas flames 

By J.-M. Samaniego 1 , F. N. Egolfopoulos 2 AND C. T. Bowman 3 


The effects of chemistry and turbulence on NO formation in oxygen- natural tur- 
bulent diffusion flames gas flames have been investigated. The chemistry of nitric 
oxides has been studied numerically in the counterflow configuration. Systematic 
calculations with the GRI 2.11 mechanism for combustion of methane and NO 
chemistry were conducted to provide a base case. It was shown that the ^simple’ 
Zeldovich mechanism accounts for more than 75% of N 2 consumption in the flame 
in a range of strain-rates varying between 10 and 1000 s~ l . The main shortcom- 
ings of this mechanism are: 1) overestimation (15%) of the NO production rate 
at low strain-rates because it does not capture the reburn due to the hydrocarbon 
chemistry, and 2) underestimation (25%) of the NO production rate at high strain- 
rates because it ignores NO production through the prompt mechanism. Reburn 
through the Zeldovich mechanism alone proves to be significant at low strain-rates. 
A one-step model based on the Zeldovich mechanism and including reburn has been 
developed. It shows good agreement with the GRI mechanism at low strain-rates 
but underestimates significantly N 2 consumption (about 50%) at high strain-rates. 
The role of turbulence has been assessed by using an existing 3-D DNS data base of 
a diffusion flame in decaying turbulence. Two PDF closure models used in practical 
industrial codes for turbulent NO formation have been tested. A simpler version 
of the global one-step chemical scheme for NO compared to that developed in this 
study was used to test the closure assumptions of the PDF models, because the data 
base could not provide all the necessary ingredients. Despite this simplification, it 
was possible to demonstrate that the current PDF models for NO overestimate sig- 
nificantly the NO production rate due to the fact that they neglect the correlations 
between the fluctuations in oxygen concentration and temperature. A single scalar 
PDF model for temperature that accounts for such correlations based on laminar 
flame considerations has been developed and showed excellent agreement with the 
values given by the DNS. 


1. Introduction 

This study is an investigation of the effects of chemistry and turbulence on nitric 
oxide formation in oxygen-natural gas flames. The choice of oxygen as the oxidizer 
is related to current interest in use of oxygen for high temperature combustion in 
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industry. In terms of NO emissions, oxygen is advantageous compared to preheated 
air in high-temperature processes such as those encountered in the glass and steel 
making industry due to the low nitrogen content. However, molecular nitrogen still 
is present in certain amounts in the natural gas (0.5% for Algerian gas to 11 %) for 
Groningue gas) and in the oxygen stream. In the latter case, the nitrogen content 
depends on the production method: from virtually 0 %t for cryogenic oxygen to 5 % 
for vacuum swing absorption. Other sources of nitrogen for oxygen-natural gas 
combustion include air leaks into the furnace. Despite the relatively low levels of 
NO emissions from oxy-combustion, ever more stringent emission standards require 
a better understanding of NO formation in oxy-flames. 

The chemical mechanisms controlling NO formation are well known (Fenimore 
1971, De Soete 1974, Miller & Bowman 1989, Drake & Blint 1991, Bozelli et al. 
1993, Bowman 1992). Nitric oxide is formed from two sources: molecular nitro- 
gen, N 2 , and fuel-bound nitrogen. In the case of oxygen- natural gas flames, the 
only source is molecular nitrogen since natural gas, whatever its origin, does not 
contain nitrogen-bound species. In this case, the three main pathways for NO for- 
mation are the Zeldovich, prompt, and jV 2 0 mechanisms, and the main pathway 
for NO destruction is the reburn mechanism. The Zeldovich mechanism is based 
on O- atom attack of jV 2 through O 4 - A r 2 — NO + A r , and is active both in the 
flame zone and the postflame zone. It is strongly dependent on temperature due 
to the high activation energy, E zeldovich ? of the jV 2 -consuming step (Ezeldovich — 
76.5 heal f mole ). A particular case of the Zeldovich mechanism is the equilibrium 
Zeldovich mechanism, where O-atoms are in partial equilibrium with molecular 0 2 . 
In this case, a global one-step reaction for NO formation can be derived with an 
overall activation energy, Eequilibrium — 138 heal j mole. The prompt mechanism 
is active only within the flame zone since it requires the presence of CH X radicals 
for consumption of molecular nitrogen, mainly through CH + iV 2 = HCN + N. 
Subsequent elementary reactions lead to the formation of N O. This mechanism is 
weakly dependent on temperature due to the low activation energy of the main iV 2 - 
consuming step (E promp t = 20 kcal/mole). The A r 2 0 mechanism is due to reaction 
of O-atoms with A r 2 in a three-body reaction, i.e. O + A r 2 -f M — A r 2 0 + M, with 
the subsequent reaction of iV 2 0 to form mainly through O + N 2 0 = NO + NO. The 
reburn mechanism is responsible for consumption of NO within the flame zone to 
produce A r 2 , and it is controlled by CH , Ci/ 2 , and CH 3 radicals. 

These various mechanisms are present in air flames and must be accounted for 
to predict accurately the level of nitric oxide emissions. The usual picture for A r O 
production in air flames can be split in two parts: 1 ) production in the flame zone 
through a balance between the prompt, A r 2 0, and reburn mechanisms - in this 
zone, the Zeldovich mechanism can be neglected; 2) production of NO in the post- 
flame zone by the Zeldovich mechanism alone. The picture is different in oxygen 
flames since in this case, the destruction rate of A r 2 is controlled by the Zeldovich 
mechanism, with a small contribution of the prompt mechanism (Samaniego et al. 
1996). The main reason is that the higher temperatures of the oxygen flame tend 
to increase the destruction rates of iV 2 , and this acceleration is more pronounced 
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for the Zeldovich mechanism due to its higher activation energy. As a result, the 
Zeldovich mechanism becomes faster than the prompt NO mechanism in the flame 
zone of an oxygen flame, while it is slower than the prompt NO mechanism in the 
flame zone of an unpreheated air flame. However, it is not clear whether the NO 
production rates can be derived from the Zeldovich mechanism alone since reburn 
may be taking place as indicated by the high HCN concentration levels observed 
in the numerical study of Samaniego et al (1996). Furthermore, it is necessary 
to clarify the role of non-equilibrium atomic oxygen since it is unclear whether the 
partial-equilibrium assumption holds (Samaniego et al. 1996). 

The effect of turbulence on NO emissions in jet diffusion flames has been studied 
extensively (Peters k Donnerhack 1981, Turns k Myrh 1991, Chen k Kollman 1992, 
Driscoll et al. 1992). The turbulent mixing process results in temporal fluctuations 
in temperature and species composition which influence the NO formation rates. 
Since the relationships between NO formation rate, temperature, and species are 
highly non-linear, NO emissions cannot be predicted from mean temperature and 
species concentrations alone. Therefore, accurate predictions of NO formation rates 
require the knowledge of temperature and composition fluctuations. Modelers have 
derived PDF formulations that account for the effect of these fluctuations on the 
turbulent NO production rate (Janicka k Kollmann 1982, Pope k Correa 1986, 
Correa k Pope 1992). The closure of the turbulent source term for NO is made 
possible through the use of a joint PDF which, in the case of Zeldovich NO with 
the partial equilibrium assumption for atomic oxygen, can be expressed as: 


= j k[0 2 ] 1/2 {N 2 }exp(-E/RT)P([0 2 },[N 2 lT)d[0 2 }d[N 2 }dT (1) 

where P([0 2 ], [N 2 ], T) is a joint PDF of oxygen concentration, [0 2 ], nitrogen concen- 
tration, [jv 2 ], and temperature, T. Usually, the nitrogen concentration is considered 
constant and [JV 2 ] is taken out of the PDF. Furthermore, in practice, [ 0 2 ] and T 
are assumed to be independent variables and the joint PDF, P, can be expressed 
as the product of two single- variable PDF’s, Pq 2 and Pt , such as: 


P([0 2 ],T) = P Oj ([0 2 ])Pt(T) 

In such a case, the turbulent NO production term becomes: 


(2) 


d[NO] 

dt 



exp {-E/RT)P T (T)dTj (3a) 


In some cases, the same term is calculated from a single PDF by assuming that 
[0 2 ] fluctuations can be neglected: 


d[NO] 


= k[N 2 )[0 2 ] 1/2 (j exp(— E / RT)P r (T)dT ^ 


dt 


(36) 
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In both cases, this final assumption has a limited domain of validity since it 
presumes that the oxygen concentration fluctuations are not correlated with the 
temperature fluctuations. This may be valid in the post-flame zone where the 
combustion process is completed, but it is highly questionable in the flame zone 
where the temperature and oxygen concentration levels are obviously correlated. 
Since the amount of NO produced in the flame zone of oxygen flames can be very 
significant, it is expected that such PDF modeling approaches are inapplicable. 

The goals of this study are to derive a simple chemical model for NO formation in 
oxygen diffusion flames and to propose a model for the closure of the turbulent NO 
production term based on a PDF approach. For this purpose, two complementary 
numerical approaches are used, one addressing the effects of chemistry, and another 
addressing the effect of turbulence. 

2. Effects of chemistry 

2.1 Numerical approach 

Concerning chemistry, NO formation mechanisms are investigated using the 
counterflow flame problem. The counterflow configuration is often used to address 
chemistry- turbulence interactions with detailed chemistry, and it is well known that 
the strain-rate modifies the chemical pathways through reduced flame temperature 
reduction and residence times (Hahn and Wendt 1981, Haworth et ai 1988, Drake 
and Blint 1989, Chelliah et ai 1990, Mauss et ai 1990, Vranos et ai 1992, Takeno 
et ai 1993, Egolfopoulos 1994a, 1994b, Nishioka et ai 1994, Samaniego et ai 
1995, Egolfopoulos Sz Campbell 1996). In such a geometry, the flow field is that of 
a strained flame where a jet of oxidizer impinges upon a jet of fuel. The oxidizer 
is a mixture of oxygen and nitrogen. Various mixtures are used: 100% O 2 and 0% 
iV 2 ; 95% O 2 and 5% iV 2 ; 21% O 2 and 79% iV 2 . The fuel is a mixture of methane 
and nitrogen. Two different mixtures are used: 100% CH\ and 0% iV 2 ; 95% CH 4 
and 5% iV 2 . 

The numerical simulation is conducted by solving the steady and unsteady equa- 
tions of mass, momentum, energy, and species concentrations along the streamline. 
Details on the set of equations and numerical method can be found in Egolfopoulos 
(1994a and 1994b), and in Egolfopoulos Sz Campbell (1996). The chemical scheme 
that is used for combustion of methane is the latest GRI 2.11 mechanism, which 
accounts for 49 species and 277 reactions. Calculations are performed without ra- 
diative losses as it has been demonstrated that radiative losses from C0 2 ? # 2 O, 
CO, and CH\ play a negligible role in counterflow oxygen flames (Samaniego et ai 
1996). 

2.2 Results 

The objective of this study is to check whether the Zeldovich mechanism alone 
is capable of predicting with sufficient precision the rate of formation of NO in an 
oxygen flame. In order to do so, two sets of calculations have been carried out, one 
with the full mechanism which includes all NO formation routes, and one with the 
three reactions of the Zeldovich mechanism, namely: 
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FIGURE 1. Computed velocity profile along the stagnation streamline of a coun- 
terflow O 2 ICH 4 flame for a strain-rate of 25 s’" 1 . 
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FIGURE 2. Computed profiles of temperature (■), oxygen mass fraction (□) and 
methane mass fraction (A) in a counterflow O 2 /CH 4 flame for a strain-rate of 25 
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0 + N 2 = NO + N 

(R i) 

N + 0 2 = NO + O 

(R 2) 

N + OH = NO + H 

(*3) 


For both mechanisms, a series of about 400 steady strained flames is calculated 
for strain-rates varying from 10 to 1000 s~ l . The strain-rate is defined as the 
maximum velocity gradient in the oxidant stream. Figures 1 and 2 depict the flame 
structure for a strain-rate of 25 s~ l , which is typical of all the strain- rates computed 
in this study. This result is obtained with the full GRI 2.11 mechanism, and it is 
almost identical for the Zeldovich calculations since the NO chemistry does not 
significantly impact the fuel chemistry. The horizontal axis is the co-ordinate along 
the stagnation streamline, with oxidant being fed at x — 0 cm and fuel being fed at 
x = 1.05 cm. The feed velocities in this case are V — 20 cm/s for both streams. The 
velocity profile exhibits a stagnation point at x — 0.74 cm (Fig. 1). The velocity 
decreases from the oxygen nozzle feed, then increases due to thermal expansion 
and decreases again until crossing the stagnation point. The velocity profile has a 
similar behavior on the fuel side, and asymmetries are due to the different molecular 
weights of 02 and CH 4 and to chemical effects. Figure 2 depicts the profiles of mass 
fractions of 02, Ci? 4 , and temperature. The maximum temperature is reached at 
x = 0.56 cm which corresponds to the flame zone. The consumption of oxygen and 
methane are almost complete at x = 0.65 cm, a location slightly on the fuel side 
relative to the maximum temperature. This can be attributed to the fact that some 
recombination reactions, which are exothermic, occur on the oxygen side. Figure 3 
depicts the profile of NO mass fraction, Y^o , in the case of the full and Zeldovich 
mechanisms, and in this case both mechanisms give similar maximum values but 
slightly different profiles in the fuel side. 

To allow comparison with well known results, a series of strained unpreheated air 
flames also has been calculated. Figure 4 depicts the maximum NO mass fractions, 
YNO,max 1 as a function of strain-rate for the oxygen and air flames, and for the 
Zeldovich and full mechanisms. In all cases, l\rvo T maz decreases with increasing 
strain- rate, and this can be attributed to reduced residence times. In the case of 
the oxygen flame, Y^o,max is well predicted by the Zeldovich mechanism alone. In 
contrast, in the case of the air flame, Yj^o.max is underpredicted by the Zeldovich 
mechanism by one to three orders of magnitude. The difference in behavior between 
the air and oxygen flames can be explained on the basis of the flame temperature 
as indicated in Samaniego et al. (1996). Oxygen flames are much hotter than air 
flames and, as a consequence, the Zeldovich mechanism, which is highly sensitive to 
temperature, is predominant over the prompt mechanism; in contrast, air flames are 
not as hot and the prompt mechanism is predominant over the Zeldovich mechanism 
at all strain-rates (Drake &; Blint 1989, Nishioka et al. 1994). 

A further comparison between full and Zeldovich mechanisms for the oxygen flame 
is performed by plotting the integral of the iV 2 consumption rate, f(— cf[iV 2 ]/cft)d:r, 
as a function of strain-rate for both calculations. The reason for looking at this 
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FIGURE 3. Computed profiles of iVO mass fraction in a counterflow O 2 /CH 4 flame 
with 5 % N 2 addition in both streams for a strain-rate of 25 $~ l . Filled symbols: 
full GRI 2.11 mechanism. Open symbols: Zeldovich mechanism alone. 
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FIGURE 4. Effect of strain-rate on maximum NO mass fraction. □: Full GRI 
2.11 mechanism for a O 2 /CH 4 counterflow diffusion flame with 5% N 2 addition in 
both streams. I: Zeldovich mechanism alone for a O 2 /CH 4 counterflow diffusion 
flame with 5% N 2 addition in both streams. A: Full GRI 2.11 mechanism for an 
unpreheated AirjCH 4 counterflow diffusion flame. ▲: Zeldovich mechanism alone 
for an unpreheated Air /C Hi counterflow diffusion flame. 
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FIGURE 5. Effect of strain-rate on integrated consumption rate of A r 2 . □: Full GRI 
2.11 mechanism for a O2/CH4 counterflow diffusion flame with 5% iV 2 addition in 
both streams. ■: Zeldovich mechanism alone for a O2/CH 4 counterflow diffusion 
flame with 5% N2 addition in both streams. 
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FIGURE 6. Profile of N 2 consumption rate in a strained O2/CH4 counterflow 
diffusion flame with 5% N2 addition in both streams for a strain-rate of 25 s~ l . 

□ : Full GRI 2.11 mechanism. ■: Zeldovich mechanism alone. 
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FIGURE 7. Profile of N 2 consumption rate in a strained OijCH\ counterflow 
diffusion flame with 5% N 2 addition in both streams for a strain-rate of 1065 s -1 . 
□ : Full GRI 2.11 mechanism. ■: Zeldovich mechanism alone. 


quantity rather than the flux of NO is that reburn, i.e. production of N 2 from 
NO and other nitrogen-containing species, can be clearly identified. Figure 5 de- 
picts that, in both cases, the consumption rate of jV 2 is positive and decreases with 
strain-rate. Furthermore, the Zeldovich approach shows two weaknesses: 1) at low 
strain-rates, it overpredicts the consumption of N 2y and 2) at high strain-rates, it 
underpredicts it. Overprediction of A 2 consumption at low strain-rates is evidence 
that reburn through the hydrocarbon chemistry is active in these cases, and under- 
prediction at high strain-rates indicates that JV 2 consumption through the prompt 
mechanism is significant. 

The discrepancies between the full and Zeldovich mechanisms may be explained 
by analyzing the profile of d[N 2 ]/dt at low and high strain-rates. Figure 6 depicts 
the profiles of d[iV 2 ]/cft for both mechanisms and the profiles of temperature, oxygen 
mass fraction, and methane mass fraction for a strain-rate of 25 s" 1 . For both mech- 
anisms, cl[N 2 ]/dt is negative on the oxygen side and positive on the fuel side. On the 
oxygen side, iV 2 is consumed to produce NO and other nitrogen-containing species 
and the two mechanisms give similar results. This indicates that A T 2 consumption 
is mainly due to the Zeldovich mechanism via the reaction O + iV 2 —► NO 4- N and 
that the prompt and reburn mechanisms are not active. On the fuel side, N 2 is pro- 
duced, which is evidence of reburn, and there is a significant discrepancy between 
the full and Zeldovich mechanisms. In this case, the Zeldovich pathway predicts 
some reburn through the reverse reaction (N + NO — * O -+* A r 2 ) but misses the 
reburn through hydrocarbon chemistry which is accounted for by the full GRI 2.11 
reaction mechanism. The picture is different at higher strain-rates. Figure 7 depicts 
the profiles of d[N 2 \/dt for both mechanisms and the profiles of temperature, oxygen 
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mass fraction and methane mass fraction for a strain-rate of 1065 . In this case, 

d[N 2 ]/dt is always negative everywhere in the flame for both the Zeldovich and full 
mechanisms, indicating that there is no reburn. However, the full GRI 2.11 scheme 
predicts higher consumption rates of N 2 everywhere in the flame zone both on the 
oxygen side and in the fuel side. On the oxygen side, most of the N 2 is consumed 
via the Zeldovich mechanism, but additional reactions are also taking place. On 
the fuel side, the full mechanism exhibits a ’’bump” that is missed by the Zeldovich 
mechanism and which is due to the ’’prompt” reaction, CH + N 2 — ► N + HCN . 

Therefore, the following picture can be drawn for oxygen flames: 1) the amount 
of N 2 consumption and, consequently, NO production decreases by about 50 % for 
strain-rates varying between 10 and 1000 s -1 ; 2) in this range of strain-rates, the 
Zeldovich mechanism predicts N 2 consumption within ± 25 % compared to the full 
GRI mechanism; 3 ) at low strain-rates, the Zeldovich mechanism overpredicts N 2 
consumption by about 15 % because it does not account for reburn through the 
hydrocarbon chemistry; 4 ) at high strain-rates, the Zeldovich mechanism under- 
predicts the N 2 consumption by about 25 % because it does not account for N 2 
consumption through the prompt mechanism. 

2.3 An overall one-step model for N 2 consumption 

The next step consists in deriving an overall one-step reaction mechanism for N 2 
consumption based on the Zeldovich mechanism alone. The proposed model (model 
1) includes reburn and considers the following reactions: 


NO + N = N 2 +0 

(fli) 

N + 0 2 = O + NO 

(R 2) 

N + OH = H + NO 

(A 3 ) 

H + 0 2 = OH + 0 

(A 4 ) 

H + OH + M = H 2 0 + M 

(A 5 ) 

O T 0 Af — 0 2 T A / 

(A6) 


The model is based on the following: 1) reaction 2 plays a negligible role (its net 
reaction rate is 1 to 2 orders of magnitude less than those of reactions 1 and 3); 2) 
atomic nitrogen is in steady-state (d[N}/dt = 0 ); 3 ) reactions 4 to 6 are assumed in 
partial equilibrium. After some algebra and considering JV-atom conservation, one 
finds: 


d[NO] 

dt 


- ~2 


d[N 2 

dt 


2 KK[h}{no] 


Kkf[N 2 }[0H}[0} 


k+{NO } + k+[OH] 


( 4 ) 


with 




[jhom 

A4 A 5 [O2] 


1 >/2 


( 5 ) 
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FIGURE 8. Effect of strain-rate on integrated consumption rate of A T 2 in a O 2 /CH 4 
counterflow diffusion flame with 5% iV 2 addition in both streams. Solid line: Full 
GRI 2.11 mechanism. • : Zeldovich mechanism alone. A: model with no reburn. 
A: model with reburn. 


[0H} = 


K 4 [H 2 0)[02] 

A’»[OJ 


1/2 


( 6 ) 


[ 0 ] 


[ 0 2 


1 !/ 2 


7) 


where bracketed quantities refer to concentrations, and where k is the specific re- 
action rate, K is the equilibrium constant, subscripts refer to reaction numbers, 
superscript ‘+’ refers to forward rate of reaction, and superscript refers to back- 
ward rate of reaction. 

An alternative model (model 2) with no reburn effects and with the partial equi- 
librium assumption for atomic oxygen can be derived by assuming k f = 0. One 
obtains: 


d[NO] 


dt 


d[N 2 


no reburn 


dt 


= 9- 


no reburn 


K, 


1/2 


[0 2 ] 1/2 [A t 2 ] 


( 8 ) 


These global one-step models are compared with the full GRI 2.11 and the Zel- 
dovich mechanisms by focusing on the evolution of the spatially-integrated A r 2 con- 
sumption rate with strain-rate (Fig. 8). This figure also depicts the values of 
f(—d[N 2 }/dt)dx of the full GRI 2.11 and Zeldovich mechanisms. As in Fig. 5, all 
models lead to a decrease of the integrated consumption rate of A T 2 with strain- 
rate. Model 1 has a behavior similar to the Zeldovich mechanism except that at 
strain-rates greater than 200 s -1 it underpredicts N* consumption. This is due 
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to the fact that for the higher strain-rates the partial equilibrium assumption for 
atomic oxygen fails. Compared to the full GRI 2.11 scheme, model 1 overpredicts 
iV 2 consumption at low strain-rates by about 15% due to the fact that it does not 
account for reburn through the hydrocarbon chemistry. At high strain-rates, it un- 
derpredicts iV 2 consumption by about 50% due to the fact that it does not account 
for the prompt mechanism and for non-equilibrium effects of atomic oxygen. Model 
2 overpredicts significantly jV 2 consumption at low strain-rates and behaves like 
model 1 at high strain-rates. The overprediction of N 2 consumption is due to the 
fact that no reburn mechanism is included, neither that of the Zeldovich mechanism 
alone, nor that of the full GRI 2.11 scheme. This shows that inclusion of the reverse 
step, N + NO — » N 2 4* O, is essential in getting correct estimates of N 2 consumption 
at low strain-rates. 


3. Effect of turbulence 

3.1 Numerical approach 

Concerning turbulence, the closure method is investigated by post-processing the 
3D data base of a turbulent diffusion flame developed by Vervisch (1992). This data 
base was the result of a direct numerical simulation of a turbulent non-premixed 
flame in decaying turbulence. The flow field was resolved accurately by solving 
the full Navier-Stokes equations. The flame was modeled by a one-step irreversible 
reaction, O + F — ► P, where O is the oxidizer, F is the fuel, and P is the product, 
and the reaction rate follows an Arrhenius formulation: 


w = kpYopYp exp (-E/RT) (9) 

where w is the reaction rate, k is the pre-exponential factor, p is the density, Yo is 
the oxidizer concentration, Yf is the fuel concentration, E is the activation energy, 
R is the gas constant, and T is the temperature. The formalism that was used 
was that of Williams (1985): the quantities were non-dimensionalized in such a 
way that the heat of reaction was expressed in terms of a temperature jump, a = 
(Tb — T u )fTb and the activation energy was expressed in terms of a Zeldovich number, 
(3 = E(Tb — T u )/RTb 2 . The values of these parameters are: a = 0.8 and (3 = 8. The 
pre-exponential factor was selected such that the initial global Damkohler number 
was equal to one. This number is defined in Chen et al (1992) as: 


Da = 


h_ 

u 0 



wdx 


(10) 


where l t is the integral scale of turbulence, ii 0 is the rms velocity, and <5// is the 
flame thickness. The Damkohler number expresses the ratio between the initial 
eddy turnover time and a characteristic chemical time. Furthermore, the molecular 
weights of the fuel and the oxidizer were taken to be equal. 

At time t = 0, the flow field was initialized with a given spectrum for turbulence 
and with a planar laminar diffusion flame located at the center of the computational 
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domain. The calculations were carried out on a 129 x 65 x 65 grid. The resulting 
data base is composed of a series of eight different time intervals: they correspond 
to times 0, 2.1, 4.4, 6.8, 9.2, 10.0, 14.7, 17.2, where time is non-dimensionalized by 
the acoustic time, t a — Ljc ( c is the speed of sound in the fresh reactants and L is 
a reference length). 

The approach for testing various modeling formulations consists of computing 
at each time-step the volume- average of d[N 2 ]/dt and comparing this value with 
that given by the various closure models. To compute d[N 2 ]/dt, a chemical model 
is needed, and we use model 2 in this section. The reason for not using model 1 
is that it is necessary to have the NO concentration for this model. However, the 
data base does not provide this information, but provides the oxidizer mass fraction 
and the temperature field, and, therefore, only model 2 can be used. Consequently, 
reburn is not accounted for. Since the objective is to assess closure formulations, it 
is expected that model 2 is sufficient. In this case, the consumption rate of NO is 
expressed from Eq. 8 as follows: 


d[NO ] 
dt 



P (-Ta/T) 


( 11 ) 


where YofT is assumed to be proportional to [0 2 \, 1 /T is assumed to be pro- 
portional to [jV 2 ], and T a is the activation temperature. In this data base, T is 
non-dimensional and varies between 2.5 and 12.5; therefore, T a is non-dimensional 
and is set equal to: 


rp E Equilibrium^^ / rnoV) 12.5 (1^) 

R(cal/mol/K) X T t (A') 

Here, we take T& = 3000 K in order to reproduce the temperature sensitivity of 
N 2 consumption occurring in an oxygen-methane flame, and this leads to T a = 
287.5. Consequently, the exact turbulent production of NO in the direct numerical 
simulation is taken to be: 


d[NO] 

dt 


(exact 


N x NyN; 


■'E Nz ,Ny,Nz^^exp(-T a /T) 


(13) 


where N XJ N y and N z are the number of points in directions x, y and 2 , respectively. 

3.2 Results 

The PDF closure models of Eq. 3a, referred to as JPDF for joint PDF, and Eq. 
3b, referred to as SPDF for single PDF, are assessed by comparing their predictions 
with the actual NO production rates as estimated from Eq. 13. The predictions 
by JPDF and SPDF are computed as follows: 


d[NO] 

dt 


JPDF 



200 


J.-M. Samaniego, F. N . Egolfopoulos & C. T. Bowman 



time 


FIGURE 9. Time-evolution of turbulent NO production term. The terms are 
non-dimensionalized by the exact value at t = 0. • : Exact term. ▲: Joint PDF. 
A: Single scalar PDF. 



d[JVO] 
dt 

SPDF 

(146) 

Figure 9 shows a comparison between the predictions of JPDF, SPDF, and the 
exact d[NO]/dt. In all cases, the production rate of NO starts increasing and then 
decreases at later times. This can be explained as follows: at initial times, the flame 
surface is wrinkled by turbulence, leading to an increase of flame surface area and, 
therefore, to an increase in NO production; at later times, the flame is strained 
by the turbulent flow field and is locally quenched, as can be seen in Chen et ai 
(1992), thereby reducing the rate of production of NO. The JPDF and SPDF 
models overestimate significantly the amount of NO production rates, with JPDF 
being slightly better than SPDF. The reason for this overestimation is that [0 2 ] 
and T are assumed to be uncorrelated. This can be seen as follows (for simplicity, 
we consider [iV 2 ] as constant): 

~ [O 2] 1 / 2 exp(T a /T) = ([0 2 ] 1/2 rexp(-r o /:r) (15) 

where ([O 2 ] 1 / 2 )* is a weighted average of [O 2 ] 1 with: 






exp(-r g /D \ 
J v exp(-T a /T)dV) 


( 16 ) 



NO formation in oxy- flames 


201 


where V is the volume on which is conducted the averaging procedure, here, the 
computational domain. Equation 16 is a weighted average which is strongly biased 
towards the high temperatures and, therefore, strongly biased to low values of [0 2 ], 
due to chemical reaction and dilatation. Consequently, we obtain: 


([o 2 ] 1/2 r < [O2F 2 


(17) 


Furthermore, we have: 

(is) 

By combining Eq. 15, 17, and 18 and assuming that [7V 2 ] inhomogeneities do not 
modify the relationships, one finds: 


jm\o 2 ]'i 2 eM-T a iT)dv < i^j^-T.mdv 


< m [0 2 ] 1/2 f exp (-TJT)dV 

Jv 


(19) 


Therefore, not accounting for the correlation between [O 2 ] and T leads necessarily 
to an overestimation of the turbulent production term of NO. In addition, the error 
is dependent on the averaging volume as is apparent in Eq. 16. In the case of Fig. 
9, the averaging volumes over which are computed mean quantities are yz planes 
as indicated in Eq. 14a and 14b, and this corresponds to a case that minimizes 
the error. Additional computations have been carried where a single averaging 
volume corresponding to the whole computational domain has been used, and the 
corresponding estimates of d[N0]/dt by the SPDF and JPDF models are 20 times 
higher. 


3.3 A model based on laminar flame structure 

To improve the predictive capabilities of PDF formulations for NO formation, a 
correlation between [O 2 ] and T is sought assuming a laminar flame structure. The 
goal is to express {O 2 ] as a function of temperature. Figure 10 shows the correlation 
between Yo and the reduced temperature 6 at time t = 0 in the data base, where 
6 is defined by: 

T - T 

(20) 

lb lu 


The relationship has two branches, one with high values of Yo corresponding to 
the oxidizer side, and one with low values of Yo corresponding to the fuel side. 
Matching the oxidizer branch leads to the following fit: 


Ko = l.- 0.9x0 


( 21 ) 


and the fuel branch is fit by: 


Y 0 = 0.05 x 9 


( 22 ) 
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FIGURE 10. Relationship between 6 and Yo across the planar laminar diffusion 
flame at time t = 0 in the data base. 



FIGURE 11. Time-evolution of turbulent NO production term. The terms are 
non-dimensionalized by the exact value at t = 0. • : Exact term, o : proposed 
single scalar PDF model. 


Combining Eq. 13, 21, and 22 leads to the following alternate single PDF model: 


d[NO ) 
dt 


model 


1 v 1 / \/l - 0.90 v^05^\ , ^ 

N x N 9 N.^ Nl ' Ny ' Nz 2y T 3 / 2 + T 3 / 2 ) exp ^ 


(23) 


Figure 11 shows the evolution of the turbulent production term of NO for the new 
single PDF model. The agreement is very good and the predictions of this new 
model are insensitive to the choice of the averaging volume. 

This test shows that this new model is promising and can be applied in practical 
codes that use single PDF formulations for NO production. In the case of an 
oxygen-natural gas flame, reburn should be accounted for and concentrations of 
O 2 • H 2 O, and NO should be included. Therefore, correlations between [O 2 ] and 
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T, [ H 2 O ] and T, and [iVO] and T should be derived from, for example, a library of 
laminar strained flames. 

In practice, this PDF model would be able to predict the amount of NO produced 
in the flame zone; however, it would not be appropriate in the post-flame zone since 
[O 2 ] and T are no longer correlated. Therefore, two models should be used, one for 
the reaction zone, similar to that developed in this paper, and one for the post-flame 
zone of the kind of the SPDF or JPDF models. 

4. Conclusion 

The effects of chemistry and turbulence on the rate of production of nitric oxide 
in oxygen-natural gas diffusion flames has been investigated. 

It has been shown that, due to the high temperatures encountered in these flames, 
the Zeldovich mechanism is dominant and accounts for more than 75% of the rate of 
production of NO over a large range of strain-rates. The limitations of predictions 
based on the Zeldovich mechanism alone have been identified: at low strain-rates, it 
overpredicts N 2 consumption, which is the source of NO , due to the fact that it does 
not account for reburn, and at high strain-rates it underpredicts N 2 consumption 
due to the fact that it does not account for the prompt mechanism. Despite these 
shortcomings, the Zeldovich mechanism is satisfactory for predicting NO formation 
rates within 25% and an overall one-step mechanism that accounts for reburn and 
which assumes partial equilibrium of atomic oxygen with O 2 has been developed. 
This simple model gives the same results as the Zeldovich mechanism at low strain- 
rates but leads to further underprediction of N 2 consumption due to the fact that 
it cannot reproduce non-equilibrium effects. 

The effect of turbulence is to generate fluctuations of temperature and species 
concentrations as well as straining effects. It has been shown by postprocessing the 
3D DNS data base of Vervisch (1992) that usual closure models based on single 
scalar PDF s or joint PDF’s lead to significant overpredictions of NO formation 
rates. The main reason is that these models do not account for the intrinsic cor- 
relation between oxygen concentration and temperature in the reaction zone of a 
diffusion flame. A new single scalar PDF model incorporating such a correlation 
has been developed. The model assumes that the relationship between oxygen and 
temperature is that of a laminar diffusion flame. The predictions of this model are 
in excellent agreement with the results from the DNS, and it is expected that it can 
be implemented in practical codes that use single scalar PDF’s for NO formation. 
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